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For a quantum computer acting on rf-dimensional systems, we analyze the computational power of circuits 
wherein stabilizer operations are perfect and we allow access to imperfect non-stabilizer states or operations. 
If the noise rate affecting the non-stabilizer resource is sufficiently high, then these states and operations can 
become simulable in the sense of the Gottesman-Knill theorem, reducing the overall power of the circuit to no 
better than classical. In this paper we find the depolarizing noise rate at which this happens, and consequently 
the most robust non-stabilizer states and non-Clifford gates. In doing so, we make use of the discrete Wigner 
function and derive facets of the so-called qudit Clifford polytope i.e. the inequalities defining the convex 
hull of all qudit Clifford gates. Our results for robust states are provably optimal. For robust gates we find a 
critical noise rate that, as dimension increases, rapidly approaches the the theoretical optimum of 100%. Some 
connections with the question of qudit magic state distillation are discussed. 

PACS numbers: 03.65.Aa, 03.67.-a, 03.67.Pp 
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I. INTRODUCTION. 

An interesting question in quantum information theory con- 
cerns the types of processes and states that must be present in 
a quantum circuit that can perform better-than-classical com- 
putation. Relevant properties that have been studied include 
entanglement (H , quantum discord \2\ and non-contextuality 
JH]. The Gottesman-Knill theorem 01 tells us that placing 
a restriction on the states and gates at our disposal can re- 
duce the power of a quantum circuit to no better than classi- 
cal. Specifically, preparation of |0) states, measurement in the 
computational basis, feed-forward of measurement results and 
unitary gates from the Clifford group (including entangling 
gates) are no more powerful than a classical computer. Collec- 
tively, we call these operations stabilizer operations. Clearly, 
the ability to create non-stabilizer states during the course of 
a computation is a prerequisite for a quantum-computational 
speed-up. 

The model of quantum computing that we adopt - perfect 
stabilizer operations plus imperfect non-stabilizer resource - 
has practical implications for implementing a quantum com- 
puter. Buhrman et al. |@| used these kinds of ideas to show 
that any qubit-based quantum computer suffering greater than 
45% depolarizing noise cannot be more powerful than a clas- 
sical computer - i.e. they proved an upper bound on the 
noise threshold for fault-tolerant universal quantum compu- 
tation (UQC). Plenio and Virmani [6] examined the effect of 
letting noise affect stabilizer operations as well as the non- 
stabilizer resource, and subsequently showed that the upper 
bound on the noise rate was on the order of 10% depending 
on the fault- tolerant scheme one intended to use. 
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The lower bound on the noise threshold for fault-tolerant 
UQC is in the region of (just under) 1% to 5% j^dl de- 
pending on whether these bounds are obtained analytically or 
via numerical estimates, and depending on particular features 
of the fault-tolerant scheme (e.g. locality constraints). The 
noise model employed for such lower bound calculations will 
also have a significant effect on the values obtained (adver- 
sarial noise models are often used in more rigorous estimates, 
for example). Typically, it seems that the bottleneck lies in 
implementing the controlled-NOT gate fault-tolerantly. In our 
model, which assumes perfect stabilizer operations and a non- 
stabilizer resource subject to depolarizing noise, no such bot- 
tleneck exists and it turns out that the lower bound coincides 
with the upper bound of 45% HHTH]. While the assumption of 
perfect stabilizer operations may seem naive, it turns out that 
many proposals for a topological quantum computer using 
anyons have exactly this property lfl4l[l5ll . In these cases, sta- 
bilizer measurements are naturally fault-tolerant, as are Clif- 
ford gates, which are implemented by braiding anyons. For 
such systems, an additional non-stabilizer resource must be 
introduced in order to achieve UQC, and since it is not topo- 
logically protected it is expected to be highly noisy. This was 
the motivation for Bravyi and Kitaev's investigation of magic 
state distillation [16]. 

For qubits, the types of resource that provably promote a 
stabilizer circuit to UQC are fairly well understood. Bravyi 
and Kitaev iflrjll . and subsequently Reichardt il3l,[l7tl showed 
that almost all single qubit non-stabilizer states could perform 
this task, when sufficiently many copies were input into a 
magic state distillation routine. The overall effect of this rou- 
tine is to output an increasingly (with more iterations) pure 
non-stabilizer state, which in turn allows for the implementa- 
tion of a unitary gate outside the Clifford group. Campbell and 
Browne have discussed the structure of distillation protocols 
lfl8ll and the impossibility of distilling some non-stabilizer 
states using known techniques |19|. Campbell I20I1 recently 
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proved the existence of catalysis-like effects for magic state 
distillation. In a previous work l2ltl we showed that, under 
depolarizing noise, every noisy non-Clifford gate either prov- 
ably enables UQC or is implementable using Clifford gates 
only i.e. there is a tight noise threshold. In the qubit case, the 
states and gates that were most robust to depolarizing noise 
turned out to be relevant for magic state distillation routines 
and this is part of the motivation for the current work - little 
has been done to date on the question of magic state distilla- 
tion for qudit-based systems. One of our results shows that, 
even for a hypothetical optimal magic state distillation routine 
for qudits, the maximum depolarizing noise rate for which a 
cZ-dimensional state could possibly be distilled is -jtj ■ 

The relevance of the discrete Wigner function (DWF) to our 
current study arises from a particular definition of DWF due 
to Galvao l22ll . which has the property that the states having a 
negative DWF are exactly those states that are non-stabilizer 
112311 . To be more precise, Galvao's construction is based on 
the DWF definition of Wooters JH and Gibbons et al. JH], 
with the additional feature that all possible legitimate DWF 
in the Gibbons et al. construction are considered simultane- 
ously. This property means that Galvao's definition contains 
redundant information in a tomographic sense (i.e. for de- 
scription of quantum states) but makes it ideal for our present 
purposes. While the interpretation of the DWF in terms of 
discrete points in phase space is a rich subject (many different 
definitions for DWF exist e.g. that of Gross i24ll ), we will 
not discuss this aspect - rather we will use this construction 
to aid in characterizing non-stabilizer states and non-Clifford 
operations. It is worth mentioning, however, that two aspects 
of non-classicality coincide - (i) the existence of negativity 
in the quasiprobability description of a quantum state and (ii) 
the existence of states that are non-classical in the sense of 
simulability via the Gottesman-Knill theorem. We prove a re- 
sult conjectured by Wootters (as noted by |27|]) concerning 
the achievable negativity of quantum states in odd prime di- 
mension. This result involves the explicit construction of so- 
called phase point operators with a certain spectrum; the more 
general problem of classifying these phase point operators by 
spectrum was recently investigated by Appleby et al. JHI. 

We begin by re-stating some relevant facts concerning qu- 
dit stabilizer operations and defining our notation. Next, we 
introduce Galvao's DWF and rewrite it slightly in terms of sta- 
bilizer projectors. The final preliminary section discusses how 
one might go about constructing a polytope whose vertices are 
Clifford gates. The table at the beginning of section [III] (Re- 
sults) represents a summary of the most robust states and gates 
that we have found. The rest of section|ni]details our method 
for ascertaining maximally robust states, facets of the Clifford 
polytope and robust non-Clifford gates respectively. 



H. PRELIMINARIES AND NOTATION. 

A. Qudit Pauli Group and Stabilizer States. 

Throughout, we alwa ys a ssume the dimension of to be a 
prime number ( 12911 and 13011 are good references for this sec- 



tion). 

X\j) = \j+lmodd) Z\j)=& j \j) (1) 

where co = exp(27ti/cQ is a primitive dth root of unity such 
thatXZ = ®~ l ZX. Further, we define 

P m =X a Z» (d>2) (2) 

P (a \ b) =i ab ^ z (d = 2) (3) 

Using so-called symplectic notation, the general form for 
multi-particle stabilizer operators with vectors x = (xi,X2, ■ ■ ■ ) 
and z = (zi,Z2, • • • ) with x i( zt £ Z d is 

P {Al) = (X x < ®X X1 -...) (Z Z| ®Z Z2 ...). (4) 

Two operators commute if and only if the symplectic inner 
product between their vector descriptions vanishes i.e. 

P{x\z) p W) = p (x'\z') p (x\ z ) if and only if 

^Xjz'i — x\zi = x- z 1 — x' ■ z = mod d (5) 

i 

The Pauli group over « qudits is given by 

g n = {<d c P {xlz) \x,zeZ n d ,ceZ d } (6) 

An n-qudit stabilizer state \\\t) is the simultaneous eigenvector, 
with eigenvalue 1, of a subgroup, Q s , of Q„. The elements of 
Q s must be mutually commuting, and Q s can be generated by 
n non-trivial elements {gi,g2, • • -gn} G Gn i.e. 

Qs = (g\,g2,---gn), \G s \=d" 

IvXvl = 4 L s 

a geg s 

where (•) denotes a generating set. 

The elements of Q n have eigenvalues from { co A " | A: e Zj} (re- 
place co with — 1 in the rest of this section for the d = 2 case). 
Measurement of an arbitrary operator A from Q„ such that the 
result is co^ is described by the following projection operator 

n [jt] = - ^/ + co^A + co" 2A A 2 + --- + co" (£/ " 1) ' : A^ 1 ) (7) 

We will have occasion to describe both single- and two- 
qudit projection operators in detail. The (O k eigenspace of a 
single-qudit operator P( a \b) corresponds to the projector 

n (fl | 6)W = -\ (i+®~%\b) + ■■■ + (o- {d - 1)k (p (a]b) ) d - 1 ) 

(8) 

which is clearly a qudit stabilizer state (there are a total of 
d(d + l) distinct single-qudit stabilizer states). Measuring a 
two-qudit Pauli operator corresponds to projecting with 

rWzki^M*] = 2 { I+(srkp (.^2\zuzi) + ■■■ 

... + ^ d -^\ p ^M^)) d ~ l ) (9) 
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This eigenspace has dimension d and can be viewed as the 
codespace of a 2-qudit stabilizer code that encodes one qudit. 

Explicit expressions for powers of a Pauli operator can be 
derived via |30| 



(A 



x\z)> 



{mx\mz 



(10) 



B. Discrete Wigner Function. 



We use the discrete Wigner function originally defined by 
Galvao which is constructed using stabilizer MUB vec- 
tors, and which was subsequently shown to have the following 
property 



Theorem 1 (Cormick et al. 112311 ) A quantum state (pure or 
mixed) has a negative Wigner function if and only if the state 
lies outside the convex hull of stabilizer states. 

The construction of this DWF relies on associating lines in 
phase space with vectors (pure stabilizer states) from a set of 
mutually unbiased bases (MUB). More precisely, for a set of 
MUB comprising d + 1 bases, each with d orthonormal states 
we have a total of d(d + 1) states labeled as |<|n) e.g. 



The mutual unbiasedness can then be summarized as 
1 

d [ 



The Wigner function for a state p at a point in phase space 
a is given by the expectation value of the so called phase point 
operator, A a , at that point 

W a = Tr(pA„) 

The operator A a is constructed using the aforementioned 
MUB vectors |(j)y) (we explain how in detail shortly), and can 
be viewed as a witness for a state being outside the hull of 
stabilizer states i.e. a state p such that Tr(pA a ) < 0, for any 
well-defined A a , cannot be decomposed as 



where the \\\f) ) are cZ-dimensional stabilizer states. The 
complete set of phase point operators A a describe bound- 
ing inequalities, defining the polytope whose vertices are the 
did + 1) stabilizer states that exist in dimension d. 

There are c/(' r/ + 1 ) distinct phase point operators, each 



of which can be associated with a vector u € Z, 
u = (hi ,U2, ■ ■ ■ ) and U[ G Zj, 



(d+l) 



i.e. 



A(u) 



'd+l 



v A-=l 



(ID 



We make the construction of these phase point operators 
more explicit in the next section. It should be apparent from 
the above construction that Tr(A) = ^ always. 



1. Restatement of DWF. 

If we desire to specify phase point operators A (u) then we 
must fix a definition for the MUB that we are using. Our def- 
inition is not much different to 112811 except for relabeling and 
different notation. It is well known that, in prime dimensions 
d, the eigenvectors of the following set of operators consti- 
tutes a MUB 

{z,x,xz,...xz d - 1 } 

Using the symplectic vector notation for stabilizer operators, 
we can identify these d(d +1) eigenvectors with projectors 
indexed by elements m; £ Z^ e.g 

{ n (o|i)[, (1 ],n (1 | 0)H ,n (1 | 1)N ,...n (1 | i/ _ 1)[M;;+l] } (12) 

Now Eq. (fTTb can be rewritten as 



Definition 1 A phase point operator indexed by u gZ^ 
corresponds to 



(d+l) 



A (") = ^ ( n (0|i)[ Ul ] 



d+l 



E n (l|>-2)[» ; ] 



(13) 



j=2 



where n^i^m are single-qudit stabilizer projection opera- 
tors. 

It is useful to have an expression for the expectation value 
of A(u) with respect to an arbitrary state p. First define the 
Pauli basis coefficients c^) for the state p 



C(a\b) = Tr ( P (a\b)P) 



(14) 



These d 2 coefficients are enough to specify an arbitrary state 
p, since the Pauli operators form a complete orthogonal uni- 
tary basis. 

Using Eq. © and Eq. ( fTOb . one can show that 



d-l 



Tr(pn (flWW ) = -£u3-**(a-2 



9k m -iq(q+l)ab 



q=0 



\(d-q)a\(d-q)b) (15) 



Substituting the projectors that are relevant to A{u) we get 



d-l 



Tr(pn (0 | 1)[ui] ) = - £<» 9"ic ( o[ (< 



q=o 

d-l 



Tr(pn (1 | 0)[H2] ) = -£a)-^c ((li _, )a 



0) 



q=0 



d-l 



Tr(pn (1|j _ 2)N ) = -Eco-^co 



-q u j a ~jq{q+i)(j- 2 ) 



(16) 
(17) 

C((d-q)\(d-q)(j-2)) 

(18) 



The last equation ( fT8l represents d — 1 distinct terms, cor- 
responding to the remaining projectors in our definition 
Eq. (fT3l , indexed by the values 3 < j < d + 1. The expectation 
value Tr(Ap) is proportional to the sum of equations ([Tol l - 
< n~8b and — Tr(pl). Each of Eq. ( TToT ) - ( TT~8b comprises powers 
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of co multiplying Pauli coefficients, c/ a \bh an( l so we can write 
Tr(Ap) in a more intuitive form involving d x d matrices, i.e. 



CO CO" 



Tr(Ap) = 



1 a U2 co"3 



"(0|0) c (0|l) 
"(1[0) c (l|l) 



(19) 



where the dot product operation between two matrices M and 
N is to be interpreted as 

d 

v / .v Y. M U N U 
y=i 

Obviously the left hand matrix of Eq. ( fT9b is a function of 
u € Moreover, its first row and column are functions 

of only mi and ui respectively, as can be seen by referring to 
Eq. ([Tol l and Eq. (fTTb . The remaining (d — 1 ) 2 elements of this 
matrix are functions of (1*3, M4, . . . , u^+i) and can be calculated 
using Eq. ( flST l. 

As a concrete example of the preceding, the expectation 
value of A(mi,W2,M3,K4) with respect to a qutrit state p (whose 
Pauli decomposition is C( a iM) is given by 



Tr(Ap) = - 2 



CO 


CO" 1 


CO 2 " 1 


CO" 2 


CO" 3 


CO" 4 


CO 2 " 2 


C0 2 " 4+1 


(j)2i<3+2 



c (0|0) c (0|l) c (0|2) 
c (l|0) c (l|l) c (l|2) 
c (2|0) C( 2 |i) C( 2 [2) 



The equivalent expression for a qubit system p and phase 
point operator A {u\, 112,113) is 



Tr(Ap) 



j_f (-1)° (-1)". 

2 2 I (_1)«2 (_1)»3 



c (0|0) c (0|l) 
C(1[0) c (l|l) 

which is probably more readily interpreted when written as 

Tr(Ap) = I (1 + (-l)"ijc+ (-l)" 2 y + (-l)" 3 z) (20) 

where x = Tr(o. v p) etc. 

The interior of the qubit stabilizer octahedron (convex hull of 
6 stabilizer states) is typically parameterized as 

W + bl + kl<i 

so we can immediately see that Tt(A(m)p) > (V u € Z,) 
defines the same region. 



2. Robustness of States to Depolarizing Noise 

If a state p is outside the convex hull of stabilizer states then 
we must have Tr(pA(t<)) < for at least one of the phase point 
operators. We define the negativity of a state as 

Definition 2 Negativity of a state, p is denoted \N(p)\: 

W(P)\ = < (21) 







Tr(pA(K))>0,VweZ; 



Definition 3 Robustness to depolarizing noise of a state p is 
denoted p*(p) : 



p*(p) = min(p) such that 



1 



(l-p)9+P 1 = Z9t\W)(W\ ^0<9,<l, L* = iJ 

where the \\f^) are d- dimensional stabilizer states. 

Lemma 2 Negativity of a state p and its robustness to noise 
are related by 



P*(P) = 1 



1 



d 2 \N(p)\ + l 



First note that since the trace of a phase point operator is 
always 4, then the maximally mixed state | satisfies 



Tr A(m) 



-p, Vu€Z' d 



d+l 



Let A be the phase point operator that minimizes Tr(Ap) for a 
particular state p, then 



Tr A 



(l-p)p+p- 



>0 
P 



-(l-p)\N(p)\ + ^>0 



=> P>1- 
P*(P) = 1- 



1 



rf 2 |^V(p)| + l 
1 

^ 2 |A?(p)| + l 



(22) 
(23) 
(24) 

(25) 



Later we will use this expression to prove the existence of 
states that can survive depolarizing rates of up to p = -j^ 
before becoming stabilizer states. 



C. Noisy Quantum Operations and the Clifford Poly tope 



We saw in section HlBl how the set of states inside the con- 
vex hull of stabilizer states can be described by bounding in- 
equalities, where each such inequality corresponds to a Her- 
mitean operator, A. We have previously argued that access to 
states from inside this region does not improve the computa- 
tional power of a quantum computer that can only implement 
stabilizer operations (albeit perfectly). For the same reason, 
access to operations inside the convex hull of Clifford gates is 
equally unhelpful. The structure of the qudit Clifford polytope 
is of inherent interest, but we are especially interested in what 
it tells us about operations that are highly robust to noise. We 
denote a noisy (depolarized) version of non-Clifford gate, U, 
as £ u i.e. 



£ v :p 



where p' = (1 - p)UpW +p- 



(26) 
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Similarly to how we defined states' robustness to noise in the 
previous section, we define the robustness of a non-Clifford 
gate U as the minimum noise rate that takes it inside the Clif- 
ford polytope 

Definition 4 Robustness to depolarizing noise of an operation 
U is denoted p*(U) : 

p*{U) = min(p) such that 

(i-pWpu'+p^Y^qiQpcJ (o < q , < 1, £$, = i\ 



where the C; are Clifford gates. In other words, p*(U) is the 
noise rate at which E\j enters the Clifford polytope. 

Finding the bounding inequalities for the convex hull of a 
set of vertices is known as the halfspace enumeration problem. 
The reverse question, finding the vertices given a set of in- 
equalities is known as vertex enumeration. Both the halfspace 
and vertex enumeration problems get rapidly more difficult as 
the size of the problem increases e.g. computing the convex 
hull of n points in if -dimensional space requires 0(nW ) time 
1I3II1 . Various software implementations exist to solve these 
kinds of problems (e.g l32l l33tl ). In order to find the con- 
vex hull of Clifford gates using this software, one must first 
describe the gates using an appropriate real, linear representa- 
tion. For example, it is possible to represent each Clifford gate 
C € SU(d) as an orthogonal rotation matrix Rc G SO(d 2 — 1 ) 



13411 . Alternatively, one can represent the Choi-Jamiolkowski 
state corresponding to each Clifford gate by using the well 
known construction for power-of-prime MUBs (e.g. ll36lo . 

Using one of the aforementioned software packages, 
Buhrman et al. @| derived a complete list of facets (non- 
redundant halfspace inequalities) for the polytope whose ver- 
tices are the 24 single-qubit Clifford gates - the so-called 
Clifford polytope. Doing so allowed them to find the non- 
Clifford gate that was most resistant to noise. One of the 
main goals of this work is to examine the structure of the 
qudit version of Buhrman et al.'s Clifford polytope |2|, i-e. 
the structure in (d 2 — l) 2 -dimensional space whose vertices 
comprise the set of c/ 3 (d 2 — 1 ) single-particle Clifford oper- 
ations on if -dimensional systems. Unfortunately, because of 
the aforementioned computational complexity of the halfs- 
pace enumeration problem we were unable to complete the 
calculation that would provide all the facets of the qudit Clif- 
ford polytope for any d > 2. However, building on our pre- 
vious work, we predicted that the facets would have a cer- 
tain structure and could be logically deduced. Checking that 
a conjectured facet is indeed a true facet is not computation- 
ally difficult, and so we were able to verify our conjectured 
facets, and consequently find a large set of distinct facets for 
d = 3,5,7. There is precedence for deriving facets of a poly- 
tope in quantum information theory, typically in the setting of 
non-locality, whereby the vertices correspond to measurement 
outcomes of different parties and facets are tight Bell inequal- 
ities (see e.g. I35TD . 



1. Checking Conjectured Facets 

A polytope V, in D dimensions, can be uniquely charac- 
terized by a finite set of bounding inequalities called facets, 
{Fj,fi}, where f) £ M. D and f G R. We can test membership 
in V of some vector X by taking the dot product with these 
facet-defining vectors Fj 



XeP 



X Fi< f V/ 



Let us denote the vertices of this polytope as Vj - again there 
are a finite number of these. 

The Clifford polytopes that we study will live in a space 
of dimension D = (d 2 — l) 2 , where d is the dimension of the 
qudit system that is under consideration. If we construct a 
conj ectured facet {G,g}, then for it to be a true facet (see e.g. 



conject 
JUS) it 



must satisfy 



Condition 1: 



Wj, VjG<, 



and, defining 



{Vj} = {Vj\Vj-G = g} 



we must also have 



Condition 2: rank 



<- v 2 

V 



= D 



i.e. we must have D linearly independent vertices contained 
in G. Later we will rewrite this second condition simply as 
rank{y,} = D. 



III. RESULTS 

A. Summary of Results 

Table I summarizes the robustness of non-stabilizer states 
and non-Clifford gates to depolarizing noise. The rest of this 
results section addresses how these values were determined. 
We first deal with the left half of the table (robust states) be- 
fore moving on to the right half (robust gates). 



B. Robust Qudit States. 

This subsection comprises the proof of the following theo- 
rem and a discussion of some of its implications. 

Theorem 3 For all odd prime dimensions, d, there is a family 
of pure states, {|v)}, which require a depolarizing rate of p = 



d+l 



in 



(l-p)\v){y\+p- 

to become a stabilizer state. Furthermore, these states are 
maximally robust amongst all states in dimension d. 
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5 ~ 8^% 
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? 


d = 7 


I = 87.5% 
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97.6% 


? 


d 


d 
d+T 


/ 







TABLE I. Summary of results: The maximal robustness to depo- 
larizing noise of non-stabilizer states p*(p), and robustness of non- 
Clifford gates p*(U). 



In other words Tr(Ap) is minimized by finding the mini- 
mal eigenvalue of A and using the corresponding normalized 
eigenvector for p (this is actually a common technique to find 
a state that minimizes or maximizes an operator). To complete 
the proof that |A^(p)| = X is achievable, we need to show that 
there is always a phase point operator, in odd prime dimension 
d, with at least one eigenvalue A-i = — X. 

Bandyopadhyay et al. |[36ll showed that the eigenvector 
(with eigenvalue Or) of XZ b is given by 



v d m=0 

d-l 

where s m = V q 

q=m 



CO b Y m \m) 



(36) 
(37) 



We will find a family of states for all prime dimension d 
that are maximally robust to depolarizing noise. Using the 
previously-derived relationship between robustness (/?*) and 
negativity (\N\), 



P*(P) = 



1 



d 2 \N(p)\ + l 



(27) 



this amounts to finding states that have maximal negativity in 
terms of their discrete Wigner function. It was conjectured by 
Wootters that the maximal negativity achievable by quantum 
states in odd prime dimension is |A^(p)| = X. Recall that the 
Wigner function of a state is its expectation value with respect 
to a phase point operator A (u). We can use this to show that 
negativity is bounded above by |iV(p)| < ^ 



W„(p)=Tr(A(«)p) 



(28) 



if d+l N 

- I Tr(pn (0 | 1)[ui] )+ £Tr(pn (lU _2) W )-Tr(pI[) 



1 

d 



7=2 

X>-i) (0 < «i < 1) 



min [W u {p)} > -- 
a. p d 

|iV(p)| < i 



(29) 
(30) 

(31) 
(32) 



Next, we show that this limit is achievable in all odd prime 
dimensions. 

Casaccino, Galvao and Severini ll27ll proved that the quan- 
tity 



min [Tr(Ap)] 



(33) 



is minimized by setting 



p = |vi)(vi| where A = £ X k \v k )(v k \ (34) 

k=l 

and Xi < Xz ■ ■ ■ < Xd (35) 



Since 



l¥o)(¥ol - n (i|fc)[o] 
and \k){k\ =n (0 |i )[fc] 



(38) 
(39) 



we have the following matrix form for the phase point opera- 
tor corresponding to u* = ( , 0, ... 0) 



A(u*) 



mm 



^Ll¥o)(^|-I 
" \ " b=o ) 

\{\mm\A I C0^»)|m)(„|-I N 

" \ a b,m,n=0 J 



d-l 



d-l 



where s m = V q and s„ = ^ q 

q=tn q=n 

For a moment let us concentrate on the matrix 



d-l 



d-l I y d-l 



- £ co^»- s '»>|m)<«|= £ -£co 

" b.m,n=0 m,n=0 \ " b=0 



m) (n | 



(40) 



Use the facts that 



d-l 



CO 



b{s„-s m ) 



6=0 



: s„ ^ s„ 

1 . S fl Sj, 



d-l l 

and s m = ^ q = — (m — m~) mod d 



q=m 



to show that non-zero coefficients occur only when 
m{m — 1) = n(n — 1) i.e. when 

n = m mod d 
or n = 1 — m mod d 

Taking all of the above into consideration, the matrix 
TiY^lm n=Q^ b ^ S " ~ Sm '\>n){n\ takes a particularly simple form - 
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it is block diagonal with a 2 x 2 submatrix of ones and an- 
other submatrix with ones on the diagonal and antidiagonal. 
For example for d = 7 we have 



d-\ 

L > 

b,m,n=0 



b(s n —Sm) 



m) (n\ 



( l 1 

1 1 


0\ 






1 1 





10 10 





10 





10 10 


\0 


1 1 / 



(d = 7) 



Note that the diagonal and anti-diagonal of the second subma- 
trix intersect at m = 

Returning to the phase point operator A(m*) we see that sub- 
tracting the identity and adding in \ ^-){ c ^- \ gives us a ma- 
trix of the form 



A{u) = - 



/0 1 
1 


0\ 






1 





1 





10 





10 


\0 


1 0/ 



(d = 7) 



i.e. it contains two counter-identity matrices of sizes 2x2 and 
(d — 2) x (d — 2) respectively. It is known that |_fj of the 
eigenvalues of an n x n counter-identity matrix are — 1 (and 
the rest are 1). In our case this means that + 1 = 
eigenvalues of A(w*) are — k. Each corresponding eigenvector 
|Vjt) (properly normalized) is a state that is maximally robust 
to depolarizing noise. 



A(0 = £^|v*)(v*| (ki<X2...<l d ) 

k=l 

{ d-l . 
So, for k e < 1,2... — - — ' 

Tr(AK)|v,)(v,|) = -i 
=>/»*(|v*)<v*|) = l 



d+l 



(41) 

(42) 
(43) 



It was noted in [28] that each A is part of a subset of phase 
point operators of size d 2 , all of which have the same spec- 
trum. We can generate the remaining d 2 — 1 remaining phase 
point operators related to A(w*) by applying each of the d 2 — 1 
non-trivial Pauli operators to it. Taking the eigenvectors, 
with eigenvalue — k, of all phase point operators related to 



A(m*) gives us a total of 



d 2 {d- 



1) 



maximally robust states. 



Let us denote by |vi) the following eigenvector of A(u*) 



Vl = 



|0> — 1 1) 
V2 



(44) 



We show in the appendix that this vector is an eigenvector 
of a qudit Clifford operation. The so-called |r)-type magic 
states for qubits (which are known to be distillable for less 



than about 35% depolarizing noise) are qubit states that are 
maximally robust to depolarizing noise, and are also eigen- 
vectors of Clifford operations. It is an interesting open ques- 
tion whether the maximally robust qudit states, |vi), that we 
have just described are distillable by some qudit magic state 
distillation routine. 



C. Method for Deriving Facets of the Clifford Polytope. 

We describe a method of deducing facets of the Clifford 
polytope by using a simple argument: 

If some two-qudit state x is acted on by sta- 
bilizer operations only (actually this process is 
the decoding for a two-qudit stabilizer code), and 
the output of this process is a single-qudit non- 
stabilizer state p, then the original two-qudit state 
x cannot have been a stabilizer state. 

Note that the converse is not necessarily true. 

Definition 5 A Choi-Jamiolkowski state corresponding to a 
unitary operation U on a qudit state is denoted 



|/ t ,) = (I 0C /)£^ 
j=o Vd 



(45) 



Definition 6 A weight-2 two-qudit x state is one for which 

c (xiMzizi) = l f Xl =Zi=0orx 2 =Z2=0 
where c {xxMzuZt] = ^{P\ xlMz ^) 



i.e. states X are those for which a local Pauli measurement 
must have zero expectation value. 

Since the identity coefficient c (o,o|0,0) i s identically 1, we 
have potentially (d 2 — l) 2 non zero coefficients. The corre- 
sponding weight-two Pauli operators form a unitary orthogo- 
nal basis for the space of Choi-Jamiolkowski states. 

The set of states {|/c,-)}> where {C\,C2 ■ ■■C ( pt d i_i\\ are all 
the Clifford gates for a single qudit, form the vertices of the 
Clifford polytope. The noisy operation Su (that we want to 
check for membership of the Clifford polytope) will be en- 
coded in the weight-two operator x via 



(I® ft/) 



d-\ 

; E \im\ 

j,k=0 , 



^2 °( x l >*2 \Z\ ,Z2l"{x\ |Z1 ,Z2) 

(46) 



If we somehow had access to a complete facet descrip- 
tion, J 7 , of the Clifford polytope (where the individual facets 
are written as Hermitean operators, F, in d x (/-dimensional 
Hilbert space) then we could immediately test whether or not 
£ u is a Clifford operation via 

3F such that Tr(Fx) < <*=S> £ v is non-Clifford (47) 
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In the absence of a complete facet description, we instead cre- 
ate operators for which the forward implication is necessarily 
true. In turn, we show that these operators are indeed facets - 
although it is possible our method does not enable construc- 
tion of every F € T . 

Definition 7 A Clifford witness W is a Hermitean operator 
satisfying 



Tr(Wx) < => £ u is non-Clifford 



1 



d-l 



\ ~ d £ \jj)(kk\ 



The general idea for our method of constructing Clifford 
witnesses is as follows 

1 . Project arbitrary weight-2 state x into codespace of two- 
qudit stabilizer code. This produces an encoded state p. 

2. Decode p — > p using stabilizer operations only. 

3. Derive conditions on coefficients c^ Xl ^ 2 \ Zl ^ 2 ) of x such 
fhatTr(pA(i<) < 0. 

Since decoding a stabilizer code uses only stabilizer op- 
erations, we know that any x satisfying these conditions on 
ct Xl ..v 2 1 s ! ,z 2 ) must have been non-stabilizer to begin with. Each 
of these conditions on cr Xl , A , | Z1 , Z2 ) corresponds to a hyperplane 
in (d 2 — l) 2 -dimensional space - where one side of the hy- 
perplane contains only x corresponding to non-Clifford oper- 
ations. What we want is a facet - a hyperplane that forms part 
of the boundary of the convex hull of Clifford operations - so 
we need to check that conditions 1 and 2 described in Section 
III C H are satisfied. 

There are many different codes we could choose so for sim- 
plicity we concentrate on codes defined by the projector 



n 



(0,0|1,1) 



d f\j)Q\ 

7=0 
1 



where \j)(j\ = -( P{ 



-J2 ^"(0,0|1,1) j G* 7 ^(0,010,1) 



(49) 



(50) 



which projects onto the parity-zero subspace of the two-qudit 
Hilbert space. These will eventually be used to derive what 
we have called "B-type" facets in previous work. We also use 
the codes 



d-\ _ 
n (o,o|o,i)[o] = YL \ 

7=0 

where 17X71 = ^(%o|o,i),<«^(o.o|i.o)) (52) 

d-\ 



n (o.o|i,o)[o] - 

7=0 

where |7)(7I = j2< p (0,0| 1,0) ,<^ p (0 



010.1) 



(53) 



(54) 



to derive different types of facets, which were called "A-type" 
facets and "A r -type" facets in previous work |f2lll . 



Conjugating any facet by local Clifford operations must 
produce another (not necessarily distinct) facet, so each facet, 
F, that we derive can immediately be used to generate a family 
of facets 



F={F'\F' = (Ci®C J )F(C,®C j ) J( } 

VUe{l,2,...,rfV-l)} 



(55) 



(4g) all of which will have the same spectrum (by construction). 

Theorem 4 The preceding method for deriving facets of the 
Clifford polytope gives a complete description (i.e. all possi- 
ble facets), for the qubit (d = 2) case. Furthermore, the opti- 
mal non- Clifford gate and the threshold depolarizing rate are 
given by an eigenvector and eigenvalue of the Hermitean Clif- 
ford witnesses Wb respectively. 



The proof is contained in the following example and the sub- 
sequent section on finding robust operations. 



1. Example: Deriving Qubit (d=2) Clifford Polytope Facets 
Define 



n (o.o|i ; i)[o] Tn (o.o|i,i)[o 
Tr (n(o,o|i,i)[o]' t: n(o i o|i,i) 



and decoding this to a single qubit produces 



d-\ 



(56) 



Pt^P= L (m\p\n)\m)(n\ 

dec 



m,n=0 



e (0.0|/U) 



^ ^ (-l)*c(i,i|jgt)+(9*c(l,l|*+i,*) 



(-i)*c(i,i|t^)+(-0* c (i,in+i,t) 

c (0,0|/fc,i) 



(57) 



Calculating the expectation value for a phase point operator 
with respect to the normalized single-qubit p gives e.g. 



For u = (0,0,0): A(w) 



/ i i_, 

2 4 4 



(58) 



4^4 



(51) and Tr (P^(")) 



Ltez, ( c {o.o\k,k)+ c (iA\k+i.k) + (-1) c (l,l|l,£)) 



LteZ 2 c {0,0\k,k) 



(59) 



Because we are only interested in the sign of Tr(pA(t/)) the 
denominator is unimportant, and so any state x that satisfies 



Y { c (o,o\k,k)+C(i,i\k+i,k) + (-l) k C(i,i\i,k)) <0 (60) 

kez 2 

must not correspond to a Clifford operation. Since 

C(x u x 2 \zuz 2 ) = Tr (- P (!ci^[? 1 ,z 2 ) x ) (by definition) the condition on 
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X can be expressed as the expectation of an observable i.e. 
Tr(WgX) < implies a non-Clifford operation, where 



Wb - E { P (0,0\k,k) ■ 



■(-Ifpj 



kez 2 



(i,i|i,*) 



(61) 



W B 



2 2-2/ 





2 + 2/ 2 



(62) 



We omit the details, but the same kind of analysis, using the 
encoding 



p oc n. 



0.010.1) 



iitn 



(0,0|0,1) 



(63) 



and the same phase point operator A (0, 0,0), shows us that any 
T satisfying 

c (o,o|o,o) + c (o,o|i,i) + c (i,o|o,i) +c(i ,o| i ,1) <° (64) 

must represent a non-Clifford operation. This inequality im- 
plies a witness of the form 



W A 



2 





l-i 














-1- 


l + i 














-1 - 





2 



\ 



(65) 



/ 



Finally, using the closely related encoding 

p x n(o,o[i,o)[o] , En(o,o|i,o)[o] 

we find that any x satisfying 

c (0,0|0,0) + C (0,0|1,1) + C (0,1|1,0) + C (0 : 1|U) < 

must represent a non-Clifford operation. The witness is 



(66) 



(67) 



2 


l-i 








1+/ 




















-1- 








-1 -/ 


2 



\ 



(68) 



Next, we seek to ascertain whether these witnesses are tight 
against the convex hull of Clifford operations. The following 
facts, which are easily checked, ensure that the witnesses are 



tight - i.e. they are facets of the Clifford polytope. 
For all F G {A,A T ,B) 
(J Ci \W F \J Ci )>0 V/G{1,2,..,24} 
and, defining the subset 

{\JcJ} f = {\Jq) I (Jc,\W F \J c ,)=0} 

rank{|Jq)} f =9 = (d 2 -l) 2 



(69) 



(70) 



Equations d69l and d70T ) correspond to conditions 1 and 2 de- 
scribed in Sectior fTl C 1 1 It should also be noted that 



{|/Q>} 



12, 



= 14 



(71) 



(72) 



(73) 



i.e. A^-type facets contain 12 Clifford vertices each, while 
B-type facets contain 14. 

As a final step, we generate sets of facets from these canon- 
ical representatives 

w F = {w' F \w' F = {c i ®c j )w F {c i ®c j y} 

V/,;€{1,2,...,24}, F G {A,A T ,B} 

resulting in the set of facets 

S = W a UW a tUW b 

It is straightforward to verify that S contains |«S| = 120 distinct 
facets, and, using software for vertex enumeration, that the 
polytope represented by S has all 24 Clifford operations as its 
vertices - i.e S is exactly the Clifford polytope arrived at by 
Buhrman et al. @]. 

D. Robust Qudit Operations. 

If we have a facet of the Clifford polytope, how do we find 
the operation that is maximally robust to depolarizing noise 
before it gets pushed inside the polytope? Recall that a two- 
qudit state x represents a non-Clifford operation £jj if 



Tr(xW) < where 



(74) 



d-l 



x = (l®£u)[- £ \jj)(kk\ =(i-p)\j u )(j v \ +p 



M=0 



d 2 
(75) 



so clearly (Ju\W\Ju) < 0, and the more negative this quantity, 
the greater the depolarizing rate required to make Tr(xW) = 0. 

Let U be the unitary that minimizes (Ju\W\Ju), and re-scale 
W to make it have unit trace, then 



Tr \W 



(l-p)\J v )(Ju\ 



d 2 



>0 



(76) 



-(l-p)\{Ju\W\Jv)\ + -jp >0 (Tr(W) = l) (77) 



P>1 



1 



d 2 \(Ju\W\Ju) 
p*(U) = l- 



d^JulWlJ^l + l 



(78) 
(79) 
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We can use the result of Casaccino et al. to find the state 
that minimizes the expectation value with the witness 

min [Tr(pW)l is achieved by 
P 

p = |vi)(vi | where W = £A*|Vjt)(Vjt| < X 2 . . .) (80) 



Furthermore if 



|Vi> is of the form \J V ) = (I ® U) YfjZo l 3 
for some U then we are finished. If this is not the case (and 
it usually is not) then the minimal eigenvalue of W at least 
provides an upper bound on the robustness of the optimal U 
with respect to W. 



we have re-derived the result of Buhrman et al. JH] in a dif- 
ferent way; the authors of ^ formulated the optimization as a 
quadratic program and used Karush-Kuhn-Tucker conditions 
to ensure the solution was a global optimum. 

The rest of this section is concerned with describing the 
most robust gates, U opt , that we have found in dimensions 
d = 3,5,7. These are analogous to the so-called Jt/8-gate for 
qubits, insofar as they are the non-Clifford unitaries that re- 
quire the maximum amount of noise to enter the convex hull 
of Clifford gates. Intuitively it seems reasonable that such 
gates would take a relatively simple form, as they do below. 



U opt w.r.t. W 



P*{U opt ) < 1 



argmin(/£/|W|/c/) 

E/eSU(rf) 

1 



d 2 \Xi\ + l 



(81) 
(82) 



The optimal qudit (d > 3) gates that we describe were found 
by numerical optimization, so we cannot rule out the possibil- 
ity that these gates correspond to a local minimum, although 
we find this unlikely. Furthermore, we cannot claim global 
optimality with respect to the Clifford polytope without first 
having a complete description of all its facets. 

In the case of the qubit facets we derived in Section ITlIC II 
the eigenvectors and eigenvalues immediately give us the 
globally optimal gate and its threshold noise rate. Consider 
the B-type representative facet Wg of Eq. d62l , re-scaled to 
have unit trace; its eigenvectors and eigenvalues are 



Vl 



v 2 



|vi) = (/®£/)I^§ 

7=0 v 2 



p*(U) = 1 



( l=i\ 

2 




. 73 J 
fl±i\ 

2 




V 73 J 



I 3711 

' -eT 



(83) 



(84) 



U 



4|i(l-V2)| + l 





0.453 



e"S" 



(85) 



(86) 



The W A (r)-type facets have eigenvalues ^(1 ± a/3) which im- 
plies a threshold noise rate of at most 1 — 1 / \/3 w 42%. Thus 



1. Robust Qutrit Operation. 



(3) 

The facet, Wg , used for this optimization is given explic- 



itly in Appendix B. 



(3) 

U opt ■= argmin(yf/|Wg '\J V ) 

£/£SU(3) 



(Ju opt \wP\J Uoi 



(87) 
(88) 



- (3 - V3cos(^) - 6cos(5) - 3sin(^) +2V3sin(f 



U, 



opt 



\ 












2711 


e"r 


2m 





e>~~r 
■p*(U opt ) Ps 78.63% 



(89) 



(90) 



2. Robust (d=5) Qudit Operation. 



The facet, Wb, used for this optimization is given explicitly 
in Appendix B. 



Uopt ■= argmin (Ju\W B \Ju) 

uesv(5) 

(Ju opt \W B \J Uopt ) = ~ 
( 



Uopt — 



(91) 
(92) 
























-2711 








4jci 

eT 








-4711 








2lH 
5 












(93) 



20 

P*(U opt ) = - w 95.2% 



(94) 



11 

3. Robust (d=7) Qudit Operation. ACKNOWLEDGMENTS 



The facet, Wb, used for this optimization is too large to re- 
produce but can be constructed by decoding the parity code 
Eq. ( TSUI ) and testing the qudit output state with respect to the 
phase point operator A(0, 0,0, 3, 2, 4, 2, 3). 
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U op t := argmin (Ju\W B \Ju) 



U 



pi 



E/eSU(7) 












b\Ju 


ipt ) 




-0.8411 






( 

1 






































lOlD 

e v 

















6m 

















e^~ 

















1 

















6TT1 

e^~ 














l° 


1 


















opt) 



'91.63% 



IV. SUMMARY AND OPEN QUESTIONS. 



(95) 
(96) 



(97) 



(98) 



Appendix A: Robust States as Eigenvectors of Clifford 
Operators 

Here we show that the maximally negative state 

|o)-|i) 



vi = 



V2 



(Al) 



is an eigenvector of a qudit Clifford operation for all dimen- 
sions d. The fact that every Clifford operation in odd dimen- 
sion d can be associated with a matrix F £ SL(2,Z/) in addi- 
tion to a vector % £ I? d results from the isomorphism 



C^SL(2,Z d ) Kl? d , 



(A2) 



established by Appleby 13711 . where C is the Clifford group. If 
we specify the elements of F and % as 



F = 




X = 




(A3) 



By using an appropriately-defined discrete Wigner function 
(DWF), we have found a family of qudit non-stabilizer states 
that are maximally robust to depolarizing noise; i.e., they are 
the states that are farthest outside the convex hull of all the 
cZ-dimensional stabilizer states. The threshold noise rate for 
these states takes a particularly simple form and, interestingly, 
these states are eigenvectors of qudit Clifford operators. Turn- 
ing our attention to non-Clifford unitary gates, we found some 
gates which required a high rate of depolarizing before they 
enter the convex hull of Clifford gates (the so-called Clifford 
polytope). In order to find these robust gates it was first nec- 
essary to deduce facets of the Clifford polytope, and we ex- 
plained a simple procedure that successfully produced many 
distinct facets. 

An obvious open question is whether non-stabilizer qudit 
states, like those that we have discussed, can be purified using 
only stabilizer operations i.e. qudit magic state distillation. 
Another natural extensions of our work is to allow noise to 
affect stabilizer operations too, as the authors of |6J] did for 
the qubit case. The depolarizing noise model we have used 
is ubiquitous in quantum information theory because of its 
generality and tractability. Nonetheless, other noise models 
are also worthy of investigation because of their relevance in 
fault-tolerance threshold lower bound calculations, for exam- 
ple, and so we highlight this as another interesting open ques- 
tion. 



then Appleby provides an explicit description of the unitary 
matrix Cr F ,%) in terms of these elements. Initially, the Clifford 
operations we are interested in correspond to 




(A4) 



which have matrix form 

C m = L* J \ d -MJ\ K = e^) (A5) 

7=0 ' 

Conjugating this Clifford operation with the Pauli operator 

%+i)/2 1 o)=Llj + (<* + l)/2>01 (A6) 

7=0 

we arrive at another Clifford operation C' 

C = P ((d+l)/2,0f(F,x) P ((d+l)/2,0) ( A? ) 

= £ |; + (d+ l)/2)^\(d + l)/2-j\ (A8) 

7=0 



For example, in dimension d = 7 we have 



12 



C 



/ A 

' 


2711 

e~r 



















2m 

e 7 


o 


o 


() 


o 


() 

























e~r 



















6TT.1 

e~^ 



















1 



















e~~r 













I o 





e^~ 














) 



Clearly |vi) = is an eigenvector of C' with eigenvalue 

— e^r. In general when j takes on the values in the ex- 
pression for C' we get the following matrix entries 



(A9) 



w f) =w u=mi,2) 



w 



(2) 



w 



(0,0|1,1)[0] 
(1,0|0,2)[0] 



+n 
+ n 



(0,0|1,2)[0] " 
(1,0[1,1)[1] 



n 



(1,0|2,1)[2] +n(i ) o[2,2)[2] 



(2) . 



144 



n(i i o|o,i)[o]+ 
■ri(i i o|i ! 2)[i]+ 



\{\Jc,)l 
rank { |/ Cj ) j = 62 (Not a Facet) 



(B4) 
(B5) 



(B6) 



Since 



2 • 



2 



d-1, 



K ( 2 



then ^^i 1 ^ is an eigenvector with eigenvalue X ■ 



|0)k ( ~ ) (1| 



|1)k 



CO* (for some k <EZ C /) 



(A10) 



-or 



(l) 



.4' 

(1) _ 



--W 



K=(0,0,0,0) 



Kr = ( n (0,0|l,l) 



n. 



w 



A>' 



i+n(o,o|i,2)[o]+n(o,i|i I o)[o]+ 

(0,l|l,l)[0]+ n (0,l|l,2)[0]+n(o,l|2,0)[0]+ 



n 



(0,1|2,1) 



n 



(0,1 12,2)[0] ' 



|l)/3 



(1) . 



= 144 



|K>}| 

rank j |7 C .) [ = 64 (True Facet) 



(B7) 
(B8) 



(B9) 



Appendix B: Facets of Qudit Clifford Polytope 

Here we explicitly describe the Clifford Witnesses that 
were used in the proving the results contained in the main text. 
For d = 3 we list all distinct Witnesses (up to Clifford conju- 
gation) that we were able to find and note whether they are 
true facets or not. For d = 5 we give only the polytope facet 
that was relevant to U opt listed in the main text. The d = 7 
facet that was used to derive the threshold and optimal gate is 
too large to reproduce here. 



1. Qutrit Facets and Peaks 



wj- = (n( ,o|i,i)[o] ^ 



n 



(0,0|1,2)[0] +n(0,l|l j 0)[0] 



.1 + 



w 



(2) 



n (0,l|l,l)[l] + n (0,l[l,2)[2] + n (o,l|2 ; 0)[0] + 

7 

n (0,l|2,l)[l] +n( ,l[2,2)[2] - 3 1 )/ 3 

|{l^>}| = 144 
rank { 1 Jq) j = 62 (Not a Facet) 



(BIO) 
(Bll) 



(B12) 



w p =w u= m o,o) 



w 



(1) 



(n 
n 



(0,0|1,1)[0] 
(1,0|0,2)[0] ' 



■11(0,01 1,2) [o] +n ( 



n 



(1,0|1,1) 



n, 



1,010,1) 



1,0|1,2)[0]~ 



n (l,0|2,l)[0] + n (i,o|2,2)[0] - o 11 )/ 3 



(Bl) 
(B2) 



w m =wr (o,o,o,o) 

W B = ( n (0,0|l,l)[0] 



n 

n. 



n 



n 

n, 



(0,0|1,2)[0] +n(L2|0,0)[0] 



(1,2|0,1)[0] _ r 11 (l,2|0,2)[0] 

n 



(1,2|1,1)[0]+ I1 (1,2|1,2)[0]+AI(1,2|2 1 0)[0]+ 



n 

n, 



(1,2|1,0)[0]" 



(1,2[2,1)[0] +n(i j 2|2,2)[0] 



10 



l)/3 



(B13) 
(B14) 



W 



(i) 



144 



{\JCi)\ =64 (True Facet) 



rank 



(B3) 



W, 



(i) . 



{|/q>}|=150 

=64 (True Facet) 



rank 



(B15) 



13 



w (2) ^=(0,0,1,2) 



(2) 



w, 



n(o,o|i,i)[o] ■ 


+-n( 0l o|i,2)[o] - 


hn ( 


n(l,2|0,l)[2] ' 


-t- n ( li2 |(12)[l] ~ 


hn ( 


n (l-2|l-l)[0] " 


F n (l,2|l,2)[2] " 


hn ( 


n (l-2|2,l)[l] " 


"l~n(l,2|2,2)[0] " 


10 





(2) 



= 150 



(1,2|0,0)[0]~ 
(1>2|1,0)[1]- 
1,2|2,0)[2]~ 

I)/3 



\\Jq)\ =62 (Not a Facet) 



rank 



(B16) 
(B17) 



This gives a total of 8640 distinct facets and 1080 distinct 
peaks. 



2. Qudit (d=5) Facet 



The following list of 58 5-tuples defines the facet Wb for 
d = 5 qudit Clifford polytope via 
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The witnesses that contain only 62 linearly independent 
\Jq) are not facets, they are called peaks. The number of dis- 
tinct facets (or peaks as the case may be) that can be created 
by conjugating with local Clifford operations is listed below 
where e.g. the set VV^ corresponds to all facets generated by 
W. 
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= 2420 



= (d 2 - 1 ) 2 = 576 (True Facet) 
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\W B X) \ = 1728, 



|W^ 2) |=864, |n4 3) |=5184 
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